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^ , Abstract. American options are financial instruments that can be exercised at any 

O ' time before expiration. In this paper we study the problem of pricing this kind 

of derivatives within a framework in which some of the properties — volatility and 
^ , dividend policy — of the underlaying stock can change at a random instant of time, 

^ ' but in such a way that we can forecast their final values. Under this assumption 

we can model actual market conditions because some of the most relevant facts that 
may potentially affect a firm will entail sharp predictable effects. We will analyse 
the consequences of this potential risk on perpetual American derivatives, a topic 
connected with a wide class of recurrent problems in physics: holders of American 
options must look for the fair price and the optimal exercise strategy at once, a 
■ typical question of free absorbing boundaries. We present explicit solutions to the 

, most common contract specifications and derive analytical expressions concerning the 

' mean and higher moments of the exercise time. 



PACS numbers: 89.65.Gh, 02.50.Ey, 05.40.Jc 
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CJ 

^ ' 1. Introduction 

Qh' Pricing financial derivatives is a main subject in mathematical finance with clear 
implications in physics. In 1900, five years before Einstein's classic paper, Bachelier [1] 
^ ' proposed the arithmetic Brownian motion for the dynamical evolution of stock prices 
■ with the aim of obtaining a formula for option valuation. Samuelson [2J noticed 
the structural failure of Bachelier's market model: it allowed negative values for the 
stock price, what led to undesired consequences in option prices. For correcting these 
unwanted features he introduced the geometric Brownian motion. Within his log-normal 
model, Samuelson obtained the fair price for perpetual options, although he was unable 
to find a general solution for expiring contracts. The answer to this question must 
wait until the publication of the works of Black and Scholes [3], and Merton [1]. The 
celebrated Black-Scholes formula has been broadly used by practitioners since then, 
mainly due to its unambiguous interpretation and mathematical simplicity. 

It is well established, however, that this model fails to fit some features of actual 
derivatives. In particular, there is solid evidence pointing to the necessity of relaxing the 
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assumption, present in the Black-Scholes model, that a constant volatihty parameter 
drives the stock price. Many models have been developed with the purpose of avoiding 
this restrictive condition: in Merton [5] volatility was a deterministic function of time, in 
Cox and Ross [6] was stock-dependent, Hull and White [7| proposed a model where the 
squared volatihty also follows a log-normal diffusion equation, Wiggins [8] considered 
underlying and volatility as a two-dimensional system of correlated log-normal random 
processes, in Scott [9J, and also in Stein and Stein [10], the instantaneous volatility 
follows a mean-reverting arithmetic Ornstein-Uhlenbeck process, and Heston [TT] 
introduced correlation in the preceding model, just to name a few. 

Another standard limitation of the Black-Scholes formula is that it is restricted to 
European derivatives: the option can be exercised at maturity only. However, most 
of the exchange-traded options are American: they can be exercised anytime during 
the life of the contract. Once again there is a clear connection with typical problems in 
physics: the earliest analysis of the issue of pricing American derivatives was formulated 
by McKean [1^ as a free boundary problem for the heat equation. Kim [13] provided 
an integral representation of the option price but, unfortunately, we have nowadays 
no explicit expression for the American counterpart of the Black-Scholes formula. The 
kernel of the problem is in that, in general, the optimal exercise boundary is implicitly 
defined by the integral equation that determines the price of the option. Only under 
certain circumstances closed expressions for American option prices do exist: e.g. in the 
case in which the properties of the derivative lead to a constant early exercise price, as 
in Rubinstein and Reiner [H], or when the option is perpetual, as in Kim [13], and also 
in Elliott and Chan [15], where the stock is driven by fractional Brownian motion. In 
the most general scenario analytical or numerical approximate methods must be used 
instead — see, for instance, Barone-Adesi and Whaley [16], Broadie and Detemple |17j . 
Ju [18], Broadie et al [T9|, and references therein as well. 

In this article we will generalize a market model first introduced by Herzel [20] as a 
simplified version of Naik's work. Naik [21] developed a model in which the volatility can 
take only two known values, and the market switches back and forth between them in 
a random way. Herzel let the volatility jump at most once: a suitable way for encoding 
a market that may undergo a severe change in volatility only if some forthcoming event 
takes place. Herzel formally solved the problem of pricing European options if the 
market price of volatility risk was constant. The problem was revisited in [22], where 
different risk premiums were considered and some explicit solutions were found. Here 
we will tackle the problem of pricing perpetual American options within a framework 
where volatility but also dividend rate may perform a single transition at some unknown 
instant in the future. 

The paper is structured as follows: In Section [2] we present the market model and 
its general properties. In Section [3] we introduce the concept of financial derivative, and 
explore the links between finance and physics in the context of American derivatives. 
Section H] is devoted to the subject of pricing perpetual American options: we stress 
the financial interest of these ideal derivatives and emphasize the potentials of the 
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analytical expressions found. In Section O we show some illustrative examples and 
discuss their implications. Conclusions are drawn in Section [6l and the paper ends 
with two appendices. In Appendix A we revisit the problem by following a different 
approach, and Appendix B deals with the smooth pasting condition for vanilla options. 



2. The market model 



Let us begin with the general description of our set-up. We will consider a financial 
market where the non-deterministic stock S is traded. The price evolution of this stock, 
assuming that S = S(to) at t = to, is the following: 

S{t) = S{to)+ [ fi{t')S{t')dt' + [ a{t')S{t')dW{t' -to), 

J to J to 

where W{t) is a Wiener process, a one dimensional Brownian motion with zero mean 
and variance equal to t. The drift, /i, and the volatility, a, are stochastic quantities 
whose initial values are /la and aa- After to they may simultaneously change to different 
fixed values, fib and cxfe, but such a transition can take place only once in a lifetime: 

fi{t) = flah<T + fJ'bh>T, 

o"W = (^ah<T + crbh>T- 

Throughout the text will denote the indicator function, which assigns the value 1 
to a true statement, and the value to a false statement. Note that if ^ a^, we 
will have a market model with stochastic volatility. The case in which /!„ ^ /i^ has in 
principle a more flexible interpretation, although we will concentrate our attention in 
the existence of two different (continuous in time) dividend pay-off regimes: 

6{t) = 6alt<T + Sbh>T- 

These magnitudes are stochastic because the instant r > to in which the transition 
occurs is a random variable. We will assume that r follows an exponential law: 

P{to <T< t\T{t)} = 1 - e-^(*-*»). (1) 

Note that ([1]) identifies A as the inverse of the mean value of the transition waiting time 
interval r — to, IE[r — to|^(to)] = A~^. In the previous expressions and hereafter jF(t) 
represents all the available information up to time t. 

The asset behaviour may be easily visualized when we express it in terms of returns, 
R(t; to) = log(S'(t)/S'(to)), instead of spot prices. The return will follow a drifted Wiener 
process with parameters Ca) up to time r. After that time the initial Brownian 
motion freezes and a second drifted Wiener process drives the subsequent evolution of 
the return: 

R{t; to) = Ra{mm{t, r); to) + i?b(max(t, r); r), 
with Ra,b{t; t') = aa,bWa,b{t - t') + 9a,b{t - t'), and ^^,6 = f^a,b - (yl b/2. 



Volatility and dividend risk in perpetual American options 



4 



In the most common situation the stock S and its derivatives, contracts whose price 
depends upon the value of this underlying, are the only securities affected by the actual 
value of r. When part of the risk is not directly traded in the market, the market may be 
incomplete: we will not be able to reproduce the behaviour of some assets by means of a 
replicating portfolio. In our case, if we want to hedge the market exposure of derivatives 
to volatility or dividend risk we must also include in the portfolio secondary derivatives , 
derivatives with the same underlying stock but different contract specifications. The 
immediate consequence of such constraint is that the risk premium coming from r is 
arbitrary to a certain extent, because investors can evaluate it on the basis of their own 
perceptions. We will avoid entering into the discussion of the financial consequences of 



this arbitrariness now: we leave it to Appendix A where the hedging-portfolio approach 



is taken. The most relevant point to be noted here is that we can use ([T]) in order to 
define a risk-neutral measure for S, 



[s < S{t) < s + ds\J^{to)} = { e 

SV^TT(^l{t - to) 

rt-to \e~^^ [iog(s/s(fp))-fl>,f-to)(t-to)]2 ^ 

+ / du — P 2,T-Hu,t-to)it-to) K 

Jo s^27ra2(M,t-to)(t-to) J 

where we have introduced some quantities depending on the risk-free interest rate r, the 
volatilities aa^b, and the dividend pay-offs 6ay. 

Oa,b = r - 5a,b - d'if,!'^, 

- 9aU + 9b{v-u) 
u[u,v) = and 

V 

V 

The use of this measure guarantees that 

= e-^'o('-'^)'^*'5(t) (2) 

fulfils 

E[F{t)\J^{to)]=F{to) = S{to). (3) 
Note that physical and risk- free measures coincide when 6'a,b = Oa^b- 



3. American option as a first passage problem 

Options are contracts between two parties, sold by one party to another, that give 
the buyer the right, but not the obligation, to buy (call) or sell (put) shares of the 
underlying stock at some prearranged price, the strike price K, within a certain period 
or on a specific date, the maturity or expiration time T. Sometimes K is a parameter 
but in general, depending on the contract specifications, it will be a function involving 
some other constants: 

K^ = Sit)TXols^,^>j,^. (4) 
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We will use the generic sign K as a shorthand for all the contract parameters. As a 
consequence of their privileged position, option holders will only exercise their rights if 
they obtain a net benefit. In other words, we can see options as contingent claims with 
their present value determined by the discounted value of the expected profit under our 
risk- neutral measure: 

P(to, Sito); K) = E[X{S{t*y, K)e-'^(**-*°)|^(to)], 

where X{S{t)] IK) is the pay-off function and t* is the actual exercise time. 

The notation we use in the definition of the pay-off function is not incidental. We 
will concentrate our attention on those contracts for which the pay-off is a function of 
the current value of the asset, like in the case of vanilla calls (-f-) and puts (-) where 

X^{S;K,T)=max{±{S-K),0)h<T, (5) 
and K is constant. Another typical pay-off is 

X^{S-Ko,T)=Xols>Koh<T. (6) 

We can derive this pay-off from (j4]) and ([5]): it corresponds to binary or digital options. 
For the sake of simplicity we will set Xq = 1 hereafter. Vanilla and binary options will 
be the only instances we will study in practice, although some other contracts may fit 
our requirements as well. Note that this is not the case of any exotic derivative whose 
pay-off depends on the past path of the stock, like in Asian options. Lookback options 
or knock-out options. 

When the option can be exercised at the end of the contract lifetime T only, the 
exercise time is deterministic, and the option is said to be European. If the option can 
be exercised at any time before expiration it is called American, and t* becomes an 
stochastic magnitude as well. Note that the contract is always worthless after maturity: 
the option buyer must decide under which conditions the option can be optimally 
exercised before this deadline. The decision will finally depend on the present value 
of stock price S{t) and the time to expiration, T — t. The problem is thus in essence 
a typical problem of first-passage time: we must determine at what time the process 
S{t) will touch the boundary H{t). In financial language, H{t) is named the optimal 
exercise boundary, the stock price above (below) which it is better to exercise the call 
(put) than to keep the option alive. We will define the exercise time i^(to) as the first 
time the underlying crosses the threshold given that at present time, to, the spot price 
of the asset lies in the proper side of the boundary: 

t^ito) = mm{t > to;S{t) | H{t)\S{to) ^ H{to)} . (7) 

The optimal strategy for choosing the boundary function derives from the following 
constraint: the investor must settle H{t) in such a way that the value of the alive option 
equals the pay-off of the contingent claim in the optimal exercise price, 

P^{t,S = H;K) = X^{H;K). (8) 
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The condition must be fulfiled in a smooth way as well, the smooth pasting condition |21 
ia[23]: 



dP^{t,S;K) 



dS 



dX^{S;K) 



S=H 



dS 



(9) 



S=H 



in the case in which the right-hand side of the previous expression does exist. That is, 
the option price must be continuous with continuous derivative in the asset price when 
it crosses the boundary. In conclusion, the investor must compute function H{t) at the 
same time he or she evaluates the option price. 



4. Perpetual American options 

We analyse now the problem of valuating perpetual American options, i.e. when we 
have T — to oo. It can be objected that perpetual American options have limited 
practical interest since, in general, actual derivatives expire. In this sense, one can argue 
that they represent the limiting value of a far from maturity contract, and therefore they 
may help in the pricing process if the theoretical price cannot be computed [16]. On the 
other hand, the existence of a fixed expiration time is a feature which is not shared by 
systems coming from other branches of science. The results we introduce in this section 
may be thus of interest in different fields. We will stress this interpretation later. 

The major simplification that perpetual American options bring is that the value 
of the boundary must be piece-like constant, 

H{t) = Halt<T + Hblt>T, 

given that the problem is stationary. Then, since the process S{t) is continuous, we will 
have: 

P±(to, So; K) = X^{Ha; K)E[e-^(*?-*°)l,±<J^(to)] 

+ A±(/7,;K)E[e-^(*?-*")l,±>J^(to)], (10) 

with = t^(to) and 5*0 = S{to). The problem is simpler in the case of binary options 
since there is no financial reason for holding the option alive once we are in the bonus 
region. This implies that Ha^b = Ko and therefore: 

D^{to,So;Ko) = E[e-^(*<t-*°)|^(to)]. (11) 

Here it is interesting to note that the magnitude we must compute for obtaining the price 
of the option is nothing but a typical moment-generating function of the first-passage 
time interval {t^ — to): 



,(n>0). (12) 

r=0 

In order to conserve this extra functionality in the output of the problem under analysis, 
we will keep 6a and 6b unexplicit as much as we can. In this way, we can turn risk-neutral 
results into physical results merely by replacing 6a^b with 6a,b- 
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4.1. The constant case 

Let us consider in the first place the case in which drift and volatility are constant. For 
9 = 9a,b and a = aa,b the probability density function (pdf) of can be obtained by 
invoking the reflection principle of the Brownian motion: 

P{t < t^ <t + dt\J^{to)} = 

^^,{t; to, So, H,,t)dt = . '''''' e dt, 

^27ralJt - tor 



where Xa,b = log (i?a,b/ 5*0). Under this assumption, the value of (fTOj) is well 
known [151 [21]: 



9n \^«.- 

X^iH,y,K){--^] . (13) 



Ha,.. 

Here, with b[-|jF(to)] we mean that we are using ip^^i^it; to, So, Ha,b) in the computation 
of expected values. We have also introduced constants 

Z^;. = (-^..^ ± ^^l, + 2ral,) , (14) 

a,b 

which differentiate the overall properties of calls and puts since ^ p]. Note that 
we must eventually recover (fT3l) with (/3^) for the hmiting case A ^ (A ^ 00). 
When the pay-off is ([5]) we will have a perpetual American vanilla option: 



V,Hto,So;K) = ±{H,-K) (^^^ 



The value of Hi, shall be obtained by demanding smoothness of the solution, recall ([9]), 
what implies here that 

dV,^{t,S-K) 



dS 

and therefore [16]: 



= ±1, 

S=H. 



/f. = ^A-. (15) 

For put options we have guaranteed < Hh < K because < 0. For call options it can 
be proved that > 1, an therefore we will have Hh > K. In fact, for (3^ = 1, 01 what is 
the same, for = r — al/2, the value of Hb diverges: there is no optimal boundary, the 
option is never exercised and, as a consequence, it must quote as its underlying stock, 
V^l^ {t, S{t)) = S{t). Thus, if the volatility is constant and the stock pays no dividend, 
only American puts are meaningful. For them we will have (3^ = —2r/al, and therefore 



V,'{to,So;K) 



2r \S 



at So f ^ 
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4-2. The general case 

We will now return to the most general case with the previous results in mind. Note 
that ([TU]) splits into two terms, the first one will count the realizations of the process that 
reach exercise price before the change in the dynamics, and the second one collects those 
for which the optimal boundary is hit after it. The first expectation can be reduced to: 

because 

mt^<r\H4)] = i-nr< = e-^(*^-*«). 



and the pdf of t^ is tpti'^'^ ^o, Sq, Ha), provided that t^ < r. Therefore, we can adapt the 
result in ffT3l) in order to obtain 



7^ 



where we have introduced 7^ which depends on A: 

7^ = ^ {-Oa ± Vdl + 2(r + . (16) 



The computation of the second term in (fTOj) is cumbersome. We will address a 
simpler problem in the first place, the calculation of E[e~''^'^*o ~*''''lj±^^|J?-'(r)], and after 
that we will recover the true expression thanks to the following identity: 

E[e-'-(*o"-*o)i^,^j^(t„)] = E fE[e-^(*?-*«)l,±>J^(r)] J^{to)] . (17) 



In order to compute E[e' ^|jF(r)] we must take into account that the indicator 

selects those trajectories for which either S~^{t) = maxto<t'<t S{t'), the maximum value 
of the process, or S~(t) = mintQ<(/<t S'(t'), its minimum value, has not met the optimal 
boundary at or before r: S^{t) ^ Ha- In this case the problem of finding the passage 
time renews: 

E[e-^(*^-*o)l,±>J^(r)] = e-^(^-*«)E[e-^(*"M-^)|^(r)]l^,(^)^^^, (18) 

where the definition of ^^(t) follows from ([7]). 

The expected value in the right-hand side of 0181) is straightforward — cf (IT^ : 



E[e-'-(*^W-)|^(r)] =E4e-^(*^(-)--)|^(r)] = (^^^ 

Therefore, we have to obtain the probability density function of S{t), S{t) = Soe^^'^'^°\ 
under the restriction 5'^(r) ^ Ha, in order to complete the computation of (ITTl) . 
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Note that, under the risk-neutral measure, we can express the return up to time r 
as R{T;to) = ±aaW^{T;to), with 



± 





Wa{t-to) + —{t-to 

On 



a Wiener process with a drift. In these terms both conditions 5'^(r) ^ Ha are equivalent 
to demand that the maximum of the Wiener process M^{t\ to) fulfils 



M±(r;to)= max W^it';to) < ± 

to<t'<T a a 

Thus, we have to compute the joint probability density function (f)'^{m,w,t,t') of the 
two random processes, 

P{m < M^(t; t') <m + dm, w < W^{t] t') < w + dw} 
= (p'tim, w, t, t')dmdw, 

which reads: 

0^(m,w,t,t') = 

2{2m — w) I y^,,o uj , wa 
= exp { ^ ± — 

CTr, 



2{t-t') ^aa"" lot *^r- 



Therefore 



A 



'S'o 



where we have introduced another constant depending on the parameters of the problem: 

(19) 



which cancels if no change in the market is induced by the point process r. When 
Hab = Kq, as in the case of binary options, the sum of the two expectations lead to 



E[e--(*?-*o)|^(to)] 



A + 



So 
K'o 



(20) 



It is easy to check that in the two limiting cases, A and A — oo, we recover the 
right expressions. In the case in which the pay-off function is ([5]) we will have 



V^{to,So;K) = ±{Ha-K) 



±{Hb - K) 



A 



A + £± 



Do 



Ha 



So_^^' 

Ha 



(21) 



The value of Ha must fulfil the transcendental equation that follows from the smoothness 
condition ([H]), 

9± 



Haii^ - 1) - {Hb - K) 



X 



^a + a2(7± + /3±)/2 \Hb 



H. 



- = 0, (22) 
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where Hf, was introduced in (fTB]) . We have seen above how an American call option 
quotes as its underlying if the volatility remains constant and the share pays no dividend. 
This conclusion is still valid in our framework, even though aa ^ Cb- Let us consider in 
the first place Sa ^ and 6b = 0. Here = 1, Hb ^ oo, and = 6a. This reduces (122!) 
to an algebraic equation in Ha, which takes the simple form: 

7+(l + A/5.) 



Then 



V+{to,So;Ha) 



-K > K. 



So 



1 

7+ 



1 + X/6a 

for 5*0 < Ha- If we consider now the limit 6a 



Sq 
Wa 



< So, 



we can see that we recover Ha — >■ oo 



and the option is never exercised also in this case. The opposite situation, in which 
6a = but 6b 7^ 0, is much more complex. We have evidences pointing to the conclusion 
that Ha —>■ oo for any combination of the remaining parameters — see Appendix B This 
implies that the stock must pay dividends after the change or in both periods in order 
to differ call options from their underlying with certainty. 

In the case of American puts we can still consider volatility and dividend risks 
separately: Hb ^ and it is bounded since < 0, and Ha = does not solve 
equation (l22ll for A 7^ 0. Moreover, when A = the expression for Ha trivially reduces 
to: 



Ha 



-K. 



For vanilla put options, as well as for call options with 7^ 1, the most relevant feature 
of transcendental equation fl22l) is the evident fact that it formally admits multiple 
solutions. A closer analysis shows that ( l22l) may have two real and positive solutions 
at the most. In Appendix B we explain how in such a case we can always discard one 
them eventually. We also find the right expression in a practical situation involving 
put options. The most surprising fact is that there are market conditions for which 
equation ( 122|) has no real solution. Unfortunately we have found no financial reason 
that motivates this behaviour. In the next section we will consider a sample case that 
exhibits this handicap but shows no other evident singularity. 

Finally note that the above pricing expressions assume that 7^ 7^ (3^. In order to 
obtain the right limiting values when a numerical concordance in the parameters lead 
to 7^ = we have to take into account that 7^ (3^ implies that £^ —A too. In 
such a case one can easily check that, for instance, the solution in (12T!) must be replaced 
by: 

V^{to,So;Ha,K) = 



± 



{Hb - K) 



A 



So 
H 
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5. Some examples and applications 

In this section we present some practical examples that may illustrate the consequences 
of volatility and dividend risk in the properties of perpetual American options. We will 
consider, in the first place, the pay-off ([5]) and analyse the outcomes related to the price. 




Figure 1. Option prices for a perpetual vanilla call under dividend risk. Here we 
consider the implications of a possible abrupt stoppage in the dividend payment. We 
represent the price of the option, in terms of the moneyness (Sq/K), for different values 
of A. We have used typical market values for the parameters: r — 4%, Sa = 1.75%, 
(5f, = 0% and aa = db = 25%. We plot as well the pay-off function in a solid (black) 
line. 



In figured] we can see the effect of a change in the dividend policy on call prices. In 
particular we consider that the quoted firm may suddenly stop the dividend payment to 
the shareholders. The consequence is that the value of the option increases due to the 
possibility of a stoppage in the distribution of dividends. The optimal exercise price also 
increases with respect to the undisturbed one: note that A = represents the limiting 
situation in which the change is impossible. When A grows the value of the option 
attains progressively the stock price, which correspond to the limiting case A oo, 
since here 6^ = 0. For finite values of A, however, all prices must eventually converge to 
the pay-off function: we can check this property in the inset of figure [H 

The following two examples deal with the effect of a change in the volatility on 
put prices. In figure [2] we can see how the price of the put steadily increases with the 
likelihood of a sudden growth in the volatility level. The behaviour reverses in figure [31 
where a possible change of volatility would imply a reduction in the put value. Note 
that we are assuming in this figure a market set-up which lead to the presence of a 
maximum admissible value for A, A ~ 0.509, with no clear reason — see Appendix B 



Another interesting property of option put prices (not related to A) is the following: 
If we observe the insets of figures [2] and [3] we will see that in the former case all prices, 
apart from the undisturbed one, decay with the same power law whereas the in the 
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0.6 









\ ■ 




Sn/K 



Figure 2. Option prices for a perpetual vanilla put under volatility risk. We represent 
the price of the option, in terms of the moneyness (Sq/K), for different values of A. We 
analyse here the consequences of a sudden increment in the volatility of the stock. We 
have used the following values for the parameters: r — 4%, Sa — = 1.75%, era — 10% 
and (Jb = 25%. The pay-off function is depicted in a solid (black) line. 




0.5 1 1.5 2 

Sq/k 



Figure 3. Option prices for a perpetual vanilla put under volatility risk. We represent 
the price of the option, in terms of the moneyness (Sq/K), for different values of A. 
We analyse here the consequences of a severe reduction in the volatility of the stock. 
We have used the following values for the parameters: r = 4%, Sa = St = 1.75%, 
(Ta = 40% and Cf, = 25%. The previous market conditions lead to the existence of 
threshold A sa 0.509. We also plot the theoretical limiting case A — > oo for which Ha 
should equal Hb. 



latter case the exponent changes with A. In fact the first exponent is fully determined 
by the parameter values after the change. The reason lies in the fact that in this case we 
have 1/9^1 < and since < |7~| for any value of A, the exponent dominates 
the extreme behaviour. Note that this feature was subtly present in figure [T] as well. 
We finally show how we can use the result in (l20l) to obtain the moments of the 
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exercise time thanks to (IT^ . In particular, we present the mean value of this first- 
passage time when the location of the boundary does not depend on r, as in the case 
of binary options, Ha^b = Kq. We will assume that the drift terms are both positive, 
9a,b > 0, and that the process starts below the barrier. Therefore, we are computing the 
mean lifetime of a perpetual binary call: 



E[t+-to|^(to)] 



Xo 



X 



(24) 



Note that the result is specially well fitted for the asymptotic analysis of large values 
of A. If we are interested in the case in which the likelihood of a change in the market 
conditions is very small the following approximate expression becomes more helpful: 
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We must remember that the proper way of performing the computation in f|T2l) is 
by assuming that r. A, 9a-, 0^, (Ja and cTb are free parameters. Once we have formally set 
r = and obtained (12^ we have to recall that this quantity must be evaluated by using 
the physical not the risk-neutral measure. Therefore we will replace Oa^b with 9a,b- 

In figure H] we have depicted the mean lifetime of binary calls under the following 
market conditions: 6a = 1.5%, Oh = 3% and ctq = 10%. Obviously the fact that the 
drift is bigger after the transition reduces the mean time as the likelihood of the change 
increases. 




Sq/Ko 



Figure 4. Mean exercise time of a perpetual binary call under drift risk. We represent 
the mean lifetime of the option, in terms of the moneyness (So/Kq), for different values 
of A. The values for the parameters are: 9a — 1.5%, 9b — 3% and aa = 10%. 



From flMl) it becomes evident that the drift of the process is determinant in 
computing the mean first-passage time, whereas the volatility plays a marginal role: 
in fact the value of cTb appears only through 6b- This outcome is no longer true if we 
focus our interest on the second moment of the first-passage time. For the sake of 
simplicity we will assume that 9a = Ob = 9, but with aa 7^ <Jb- We can observe how, even 
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under our previous assumptions, the second moment depends explicitly on the value of 
a a and a^. 

ni4 '''' 



to)'l^(to)] 



at 



Xo 



crt 



A 



an expression that reduces to 
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(26) 



if we are concerned about small values of A. In fact, if we confront fl25l) and fl26l) 
we realize that the first correction to the mean first-passage time is governed by the 
undisturbed second moment of the process. 




0.2 0.4 0.6 0.8 



Figure 5. Standard deviation of tlie exercise time of a perpetual binary call without 
drift risk. We represent the variance of the lifetime of the option, in terms of the 
moneyness {Sq/Kq), for different values of A. The values of the parameters were 
adjusted in order to keep the drift unchanged: 9a — db — 1-5%. 



In figure [5] we plot an example of the standard deviation of the exercise time which 
will illustrate the constant drift instance. Clearly this is only feasible if the value of the 
parameters are well tuned: fia = 2%, fib = 3.5%, (Xa = 10% and ab = 20%. The outcome 
fits our anticipation, since the possibility of a larger volatility increases the uncertainty 
about the mean exercise time. 



6. Conclusions and future work 



In this article we have considered the implications that the presence in the market of 
volatility and dividend risk has for option prices. The proposed market model allows 
random changes in both the dividend-payment rate and the volatility level of stock 
shares, but in a very specific way: only one change is feasible, and the final market 
properties are foreseeable. The model, however, is rich and realistic enough to obtain 
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valid financial results. We have focused our attention on the problem of pricing perpetual 
American options: derivatives with no expiration limit that may be exercised at any 
time. The absence of maturity in a derivative departs from actual market conditions, 
but it is a well-accepted approximation used with the purpose of casting light on the 
way of solving the whole problem. From a physical point of view, the analysis of 
perpetual options may have even bigger interest than real options because, usually, 
physical systems do not disappear after a fixed time lag. 

Within this framework, we have obtained explicit solutions for pricing derivatives 
with the most typical pay-offs: vanilla puts and calls, as well as binary options. 
Nevertheless, the applications of our development are not restricted to finance since, as 
we have pointed out, binary option prices are nothing but classical moment-generating 
functions of a first-passage time. We have shown the right way of handling these 
expressions within our set-up. We have illustrated our results with a set of practical 
examples covering the major issues we have analysed. 

The most intriguing aspect of our outcomes is the possible existence, depending 
on the actual market conditions, of a threshold that limits the probability of change in 
order to have a valid price. We have found no sound clue that may explain the financial 
meaning of this behaviour, beyond the fact that in such a case the smooth pasting 
condition cannot be satisfied. We believe that this featiuc deserves deeper analysis 
in a forthcoming work. We aim to extend these results to American derivatives with 
maturity date as well. It is probable that the solution to this problem will imply the 
use of approximate procedures and numerical techniques to some extent. 
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Appendix A. The financial approach 

In this appendix we will present a supplementary approach to the issue. Both methods 
lead to the same results but differ in the way in which they are obtained. This fact can 
made a technique more suitable than the other when discussing some properties of the 
model, like the question of the completeness of the market. This alternative approach 
is based on the idea that the fair price of an option must be equal to the value of some 
portfolio made of different securities that mimics the behaviour of the derivative and 
hedges all the risk. In fact, the market will be complete if we can construct such hedging 
portfolio for every traded asset of the market. 

The first security to be included in the portfolio is the underlying asset, which will 
reproduce changes in the option price due to the evolution of the stock price S. The 
second security in our portfolio is a zero-coupon bond, a free-risk monetary asset, with 
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a market price B that satisfies the following differential equation: 

dB = rBdt. (A.l) 

A long position in this security will provide a secure resort where to keep the benefits 
of an effective investment strategy, whereas a short position in these bonds will allow 
us to borrow money when we need it. These two securities cannot counterbalance all 
the stochastic behaviour of the price of the option however: not all the influence of 5{t) 
and a{t) on the option price may be explained through S; and B is fully deterministic. 
Therefore, we need another security that can account for this contribution to the global 
risk. Nevertheless, in the most of the cases, markets do not trade such assets, a fact that 
impel us to consider the inclusion of a secondary option in the portfolio: a derivative 
of the same nature of P{t,S;K), but with a different set of contract specifications, 
Q{t, S;K'). In particular we will focus on the case in which they share the same 
expiration date but differ on the striking price. Note that we will assume as well that 
the option price is a function of the current value of the underlying. 

Let us write down P as a mixture of u shares S, (j) units of the riskless security P, 
and -0 secondary options Q: 

P^iyS + (f)B + ijjQ. (A.2) 

The variation in the value of the portfolio, due to the market evolution of its components 
and the received dividends, fulfils the following relationship: 

dP = udS + uSSdt + (f)dB + ipdQ, (A.3) 

where we have taken into account that u, 0, ip and 5 are predictable processes, and that 
we adopt a self-financing strategy, in which there is no net cash flow entering or leaving 
the replicating portfolio. This differential change must equal the expression obtained 
after applying the rules of Ito calculus on the price of the option: 

dP = dtPdt + dsPdS + ^a'^S'^dlsPdt + APdlt>r, (A.4) 

where 

AP = Pfe(t, 5; K) - P(t, S\ K), 

and Phit, S; K) is the price of the option after the jump. The differential of an indicator 
with a random variable in its argument may seem a bizarre object. However, it is 
mathematically well defined. In fact, lt>T is a submartingale under our measure and, 
by virtue of the Doob-Meyer decomposition theorem, dlt>T can be expressed as a sum 
of two terms, dlt>T — dA + dM, an increasing adapted process A, 

A{t) = X-[nm{T - to,t - to), (A. 5) 

and a cadlag martingale M, 

M{t) = 1 - AE[r - to\J^{t)] . (A.6) 

The key point is that dlt>r is a stochastic magnitude, independent of dW, which does 
not directly contribute to the variation of the stock price dS. Therefore, it is a source of 
risk that cannot be explained in terms of the random evolution of the underlying asset. 
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dtPdt + dsPdS + -a'^S^dlsPdt + APdlt>, 



1 
2 

udS + SSdt + rcpBdt + ifjdQ. (A.7) 
Now, we can proceed with dQ in an analogous way, 

dQ = dtQdt + dsQdS + -a^S'^dlsQdt + AQdlt>r, (A.8) 
where the natural definition of AQ, 

AQ = Qbit,S;K')-Qit,S;K'), 

has been used. In order to recover a deterministic partial differential equation we must 
guarantee that terms containing the stochastic magnitudes dS and dlt>T mutually cancel 
out. Therefore we must demand that 

u = dsP- tpdsQ, 

a condition named delta hedging, and also that 
AP 

^ = Ag' 

which is usually referred as vega hedging. The previous hedging conditions reduce ( ]A.7p 
to 

1 

— ( 
2 

an expression that still involves B, which is not an inner variable of option prices P 
and Q in our set-up|j| This problem can be fixed using the definition of the portfolio 
m ([Oil and the vega hedging together, 

/ AP \ AP 



dtP + ^a^S^dlsP - SSdsP 



dtQ + ^cj'S'd'ssQ - SSdsQ , (A.9) 



The replacement of (f)B in flA.9P leads to 

dtP + ^cr^S^dlsP -rP + {r- 6)SdsP = 

^ (dtQ + Icx'S'dlsQ - rg + (r - 6)SdsQ^ . 

This formula implies the existence of an arbitrary function x, which uncouples the 
problem of finding P and Q: 

X = ^ (dtP + l^x'S^dlsP - rP + (r - 6)SdsP^ . (A.IO) 

Obviously the same formula is valid for the secondary option, merely by replacing P with 
Q. This proves that the option Q removes the remaining risk and therefore completes 
the market. 

I In fact, some authors replaces t with B as the free variable. 
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The financial interpretation of x is discussed in more deptli in [SU] and [22] • Tlie 
must known facts are four: (i) x depends on how every investor measures the volatihty 
and dividend risk, what prevents us from fixing it, (ii) prior to the change it must be 
negative defined or otherwise the market will show arbitrage opportunities, (iii) after 
the change must be zero and, finally, (iv) the choice x = —^^t<T avoids the so-called 
statistical arbitrage, the growth of the expected value of the discounted price of the 
option, P = e~''*^*^*"-'P. The violation of this condition is against the capital asset 
pricing model which states that any hedged portfolio with a zero market risk must have 
an expected return equal to the risk-free rate. It can be easily shown that this choice 
for X lead to the desired property: from ( ]A.4I) and ( lA.lOl) it follows that dP can be 
decomposed in two parts, 

dP = dsPdF + APdG, 

where F was defined in ([2]) and dG = xdt + dlt>T. On the one hand F is (and must be) 
a martingale under our measure — see ([3]). On the other hand, if we set 

xdt = -dA = -Xlt^rdt, (A. 11) 

with A defined as in (IA.5I) . we will have that G equals M — see (I A. 6 II — , and therefore 
becomes a martingale as well, what proves our previous statement. Note however that 
in general x may depend on the two independent variables t and S, on market properties 
r, S, a and r, but also on the common contract specification T. 

In fact, financial arguments — see for instance [4] — lead to the conclusion that 
the natural time variable in option pricing problems is not t but t = T — t, the time 
to maturity. This convention, for instance, turns the pay-off function into an initial 
condition in European-hke contract problems P(t = 0, S = S(T); K) = X{S(T); K) and 
removes explicit temporal dependence from perpetual derivative problems, also in the 
optimal exercise barrier of American options Ymvi^ooHit) = Ha,b- This reduces (lA.lOll 
to an ordinary differential equation: 

1 d'^P dP 

r'^'ds^ + " ^^^ds " ~ = ^^■^'^^ 

The solution if the jump has taken place does not depend on x since the equation to be 
solved is [16j : 

r'^'it + " ^'^^^ - = °- ^^-^^^ 

Here x becomes meaningless and thus we can freely set x = 0. The general solution 
of flAlS]) for Sb^O reads 

where (3^ were defined in the main text — cf (|T4l) — , and constants Ci and Co can 
be obtained by demanding the solution to fulfil boundary conditions ([8]) and (Q, if 

applicable: e.g. Gi = and Gq = l/K^'' lead to a perpetual American binary put, 

Gi = {Hf, — K)/H^'' and Cq = to a perpetual American vanilla call, and so forth. 
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The solution if the jump has not taken place yet does depend on the arbitrary function 
X- The development in the main text is consistent with flA.ll|) . i.e. x = this 
case ( ]A.12[) reduces to: 

1 d'^P dP 

iflS'jS^ + (r - 5^)8"— - (r + X)P = -AP,, (A.14) 



whose general solution is 



p = C^S^^ + ^2^^" + A 



^1 .5/3^ + ^^5/3,- 



(A.15) 



_A + £+ A + 

whenever 7+ 7^ 7" and 7^ 7^ (3^. Expressions for 7^ and have been already 
introduced in the main text, in f|T6|) and f|T9|) respectively. Once again the arbitrary 
constants in flA.lSp must be evaluated on the basis of contract specifications: for example 
C3 = £+/(A + £+)irj" , Ci = l/fsTo' and C2 = Cq = reproduce the price of a perpetual 
American binary call. 

Clearly, the previous method is well suited to discuss most of the financial 
assumptions and implications of the market model but dilutes the connections of the 
issue with other branches of knowledge. The fact that the solution of flA.12p is related 
to a first-passage time problem may be not so evident at first glance, for instance. The 
way of recovering the results corresponding to the use of the physical measure is also 
a delicate question. These and other reasons impelled us to follow in the main text a 
development which is akin to physicists, and to leave the mathematical finance approach 
to this appendix. 



Appendix B. The smooth pasting condition 

In this appendix we analyse the properties of the smooth pasting condition for vanilla 
options: formula fl2^ in the main text. Let us first define the auxiliary function ^^(r/), 

with 

The study of this function is relevant in our framework since Ha fulfils Q^{Hal Hh) = 
— remember that Hh was introduced in ( ITSl) . The properties of and 7^ impel us 
to treat separately call and put optimal boundaries. For call options we will have that 
Pab — 1 and 7^^^ > 1. We will dismiss the particular case = 1 because it was 
analysed in Section HI Moreover, the only scenario in which 7^^ = 1 corresponds to 
a very specific situation: /3+ = 1 and A = 0, the jump never takes place. Therefore, 
A~^ and are bounded and positive constants, G~^{ri) fulfils lim^^o ^'''(^) = —00 and 
lim^^oo S'^iv) = —B^i and it shows a single maximum aX t] = tjm'- 
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That means that the equation Q^{r]) = may have two solutions, one solution or no 
solution at all, depending on the value of G'^^tim)- When the parameters lead to 



for any value of A we will have two formal solutions, but we can easily reject one of 
them by analysing the properties of the boundary for large and small values of A. In 
fact, condition (IB.ip is always satisfied for small enough values of A since 

V) 

The previous statement is false if = 1, i.e. when the stock pays no dividend before 
the change. In this case we will have found no parameter choice which fulfils (IB. Ill 
whenever ^ 1. In particular, it is easy (but tedious) to show that if r > and 6a = 
then 

7^-l^% A 



it J A + r + al-f+ ' 
for all A > 0. 

Let us assume again that > 1. Under these assumptions, and depending on the 
market value of the parameters, may exist a threshold for A, A, above which no optimal 
boundary can be defined: 

it-iY^ A 



TJ ''''IT 



it J \ + r + al^t(3tl2- 
In fact, for A = A, we will have a single solution: 

Jl 

This situation is also present in the estimation of the optimal boundary for puts. 
Here we have (3~^^ < and 7~ < 0, what leads to the conclusion that A~ and 

are also bounded and positive constants, that the auxiliary function Q~{f]) fulfils 
lim^^o Q'iv) = B', lim^_^oo Q'iv) = , and that 7]^, 
^ _ l7a I 

l7a I + 1 

is a global minimum of Q^{r]). The condition for the existence of a solution for the 
optimal boundary problem now reads: 

< A + r + ag|7-||/?,-|/2 ^ ^^^2) 
and, like in the call case, it will be satisfied for small values of A because 

< oo. 





+ 1 









f\Pa\ 


+1^ 
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Therefore we may have a threshold in A depending on the market properties as welL Let 
us show this behaviour through a practical example. Consider the case (3^ = 2. Here 
the constraint (IB. 21) is equivalent to demanding (r + 6*0)17"! + r > 0. This means that 
for 6a > — T we will always find a optimal boundary, whereas for 6a < —t there will be 
solution only in the range 



Note that Ha shows the right limiting properties \mi\^QHa = K\[3^ |/(|/^a I + 1) ^'^d 



We have a plausible explanation of the absence of solution for any value of A in the 
case of calls on non-paying dividend shares: the perpetual American call behaves like 
the underlying and thus there is no optimal boundary. The financial interpretation of 
the existence of threshold A, but only under some market conditions, is elusive. 

References 

[1] Bachelier L 1964 The Random Character of Stock Market Prices ed P H Cootncr (Cambridge MA: 
MIT Press) 

[2] Samuelson P A 1965 Industrial Management Rev. 6 13 

[3] Black F and Scholes M 1973 J. Pol. Econ. 81 637 

[4] Merton R C 1973 Bell J. Econ. and Management Sci. 4 141 

[5] Merton R C 1976 J. Financial Economics 3 125 

[6] Cox J C and Ross S A 1976 J. Financial Economics 3 145 

[7] Hull J and White A 1987 J. Finance 42 281 

[8] Wiggins J B 1987 J. Financial Economics 19 351 

[9] Scott L O 1987 J. Financial and Quantitative Anal. 22 419 

[10] Stein E M and Stein J C 1991 Rev. Financial Studies 4 727 

[11] Heston S L 1993 Rev. Financial Studies 6 327 

[12] McKean H P 1965 Industrial Management Rev. 6 32 

[13] Kim I J 1990 Rev. Financial Studies 3 547 

[14] Rubinstein M and Reiner E 1991 Risk 4 75 

[15] Elliott R J and Chan L 2004 Quant. Finance 4 123 

[16] Barone-Adesi G and Whaley E 1987 J. Finance 42 301 

[17] Broadie M and Detemple J 1996 Rev. Financial Studies 9 1211 

[18] Ju N 1998 Rev. Financial Studies 11 627 

[19] Broadie M, Detemple J, Ghysels E and Torres O 2000 J. Econometrics 94 53 

[20] Herzel S 1998 Int. J. Theoretical and Applied Finance 1 487 

[21] Naik V 1993 J. Finance 48 1969 

[22] Montero M 2004 Eur. Phys. J. B 42 141 

[23] Shackleton M B and S0dal S 2005 Economics Letters 89 200 

[24] Wilmott P 1998 Derivatives (Chichester: John Wiley & Sons) 




The value of Ha, when it exits, follows: 




